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We consider a phenomenological continuum theory for an extensile, overdamped active nematic 
liquid crystal, applicable in the dense regime. Constructed from general principles, the theory is 
universal, with parameters independent of any particular microscopic realization. We show that it 
exhibits a bend instability similar to that seen in active suspensions, that leads to the proliferation of 
defects. We find three distinct nonequilibrium steady states: a defect-ordered nematic in which +| 
disclinations develop polar ordering, an undulating nematic state with no defects, and a turbulent 
defective nematic. We characterize the phenomenology of these phases and identify the relationship 
of this theoretical description to experimental realizations and other theoretical models of active 
nematics. 


Liquid crystals are anisotropic fluid mesophases 
that exhibit broken rotational symmetry and have 
been extensively studied and extremely useful for 
many years [1]. The study of topological defects in 
the orientational order in these systems, which typ¬ 
ically occur under driving, has had a central role in 
developing our understanding of the material prop¬ 
erties of these systems [2, 3]. Active liquid crys¬ 
tals are anisotropic fluid phases that are driven at 
the scale of the microscopic nematogen. These mi¬ 
croscale internal forces give rise to spontaneous de¬ 
fect nucleation, and novel defect dynamics [4-6]. 

An experimental realization of a two-dimensional 
active nematic with rich phenomenology is a sys¬ 
tem of cytoskeletal filaments driven by motor pro¬ 
teins [4-6] confined to a fluid interface. This novel 
nonequilibrium system has been shown to have tran¬ 
sient self-propelled defects and emergent defect or¬ 
dering and has triggered much theoretical effort to 
understand its dynamics [6-18]. One outcome of the 
theories is to show how fluid mediation, in the form 
of active and passive backflow, can lead to the forma¬ 
tion and propulsion of defects. However, defects also 
arise in active nematic systems in which fluid medi¬ 
ation plays little or no part, such as vibrated mono- 
layers of granular rods [19], epithelial cell monolayers 
[20], and elongated fibroblasts [21]. 

In this work we develop a phenomenological con¬ 
tinuum theory that describes the dynamics of an 
overdamped active nematic, and is applicable to all 
systems in this symmetry class in two dimensions. 
The nonequilibrium steady-states that we find are (i) 
a defect-ordered nematic that exhibits emergent po¬ 
lar ordering of +4 defects, (ii) a defect-free undulat¬ 
ing nematic, and (iii) a defective, turbulent nematic 
(see Fig. 1). Further, we identify the relationship of 
existing theories to this framework. 



FIG. 1. (color online) Heatmaps of the degree of order 
(S/So, colorbar on the upper-right) showing the nonequi¬ 
librium steady-states. Top Left Panel: The defect- 
ordered state, with a histogram inset showing the sharp 
polar ordering in the orientation of + § defects. The 
lines, showing the direction of order, highlight the ex¬ 
tended trails left by the motion of these defects. Top 
Right Panel: The undulating nematic state is highly or¬ 
dered (S ~ So) but the direction of order undulates. 
Bottom Panels: Defective states at low activity (left, 
'f) < 1) and high activity (right, ^ > 1), with 2x- 
magnified regions (lower left), and heatmaps of the vor- 
ticity (in the upper left, colorbar on the lower-right). 
Scale bars are on the top right. 


Theoretical Framework: An equilibrium ne¬ 
matic is described by the well-known Landau-de 
Gennes free energy that is a functional of the den¬ 
sity p(r , t) and the nematic order tensor, Q(r, t), as¬ 
sociated with rotational symmetry breaking [1]. Its 
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dynamics is given by gradient descent on this free 
energy landscape: ‘Model A’ dynamics for the direc¬ 
tor, d t Q(f) = —y -1 ? and ‘Model B’ dynamics 

for the density, d t p(r) = — 7 /_1 V 2 , as density is 
conserved. If this system is in an imposed flow then 
the director is convected and rotated by that flow, 
and the dynamical equation becomes 

(d t + u- V)Q = — 7 _1 ^- + Qfl-nQ 
+AE +A(QE + EQ) r 

where ft is the vorticity tensor, E is the strain- 
rate tensor associated with the flow, A is the flow- 
alignment parameter, and (A) r denotes the trace¬ 
less version of A (e.g. A — \Tr{ A)I). 

In the case of an active liquid crystal, internal 
stresses due to the forces exerted by the particles 
themselves give rise to self-generated flows. We pos¬ 
tulate that an active liquid crystal undergoes gradi¬ 
ent descent dynamics in the local rest frame of this 
self-generated flow arising from the activity. The 
active stress has the same symmetry as the nematic 
order (a oc Q) and hence the self-generated flow will 
be proportional to the force density: pu = F/£,t = 
—AcVQ. Similarly the vorticity will be proportional 
to the torque density: puo = t/£r = A#(V x VQ), 
where £’s denote friction coefficients. Using the form 
of the flow and vorticity arising from the activity, the 
dynamical equations for an active nematic take the 
form 

dtQ = —7 _1 ttt + — VQ] • VQ - ^Q(VV : Q) 

0(4 p 2p 

-^(Q-VV-Q-VyQrQ) 5 (1) 

where the coefficients Ac, Aand A^ control the 
strength of active convection, active torque, and 
flow-alignment respectively, A : B = AijBij , and S 
denotes symmetrization (e.g. (A) 5 = ^(A +A T )). 
Further, we take the dynamics of the density to be 
of the form 

dtp = DV 2 P + DqVV : Q (2) 

which incorporates the self-generated flow through 
the curvature induced density flux [22, 23], con¬ 
trolled by Dq. 

Equations (1) and (2) are the dynamical equations 
of an overdamped active nematic which will be stud¬ 
ied in this work. Before we proceed with the analysis 
of our theory, we make the following observations in 


order to place this model in the context of other the¬ 
ories in the literature of active liquid crystals: 

(1) Suppose our active fluid had polar symmetry, the 
self generated flow will be proportional to the polar 
order parameter, pu = F/^t = A P. Postulating a 
gradient descent dynamics, as above, for an active 
polar fluid results in the same dynamics as described 
in the seminal works of Toner and Tu [24, 25]. In 
this sense, this theory is a generalization of their 
approach to active nematic systems. 

(2) Existing theories of active nematics consider a 
coupled set of equations of the nematic order param¬ 
eter and the active induced flow which arises from a 
Stokes equation (such as [7-13, 18]). Suppose there 
exists a screening mechanism such as confinement 
to 2D, then eliminating the flow field in terms of the 
active stress yields Eqs (1) and (2) above, but with 
Ac = A r i.e., a Galilean invariant version of our the¬ 
ory [18, 26-28]. This is an important distinction of 
our work. In over damped systems Galilean invari¬ 
ance is broken by the medium through which our 
particles move (as first pointed out in [24]) and this 
is accounted for in our theory. 

(3) Theoretical work on active nematics was pio¬ 
neered in [22, 23], and have subsequently been shown 
to have giant number fluctuations, phase separation, 
and band formation near the critical density [19, 29- 
40]. Our Eqs (1) and (2) reduce to this description 
when A# = Ac = A^ = 0. Hence our theory can be 
considered a generalization of these previous works. 

Parameters of the theory: The equilibrium 
contributions to the dynamics of the order parame¬ 
ter take the form 

[dtQij\Eq = D r [a — f3TrQ 2 }Qij + 2 DE^7 2 Qij 
- {[diQki]djQki) r ) 
+D p (d i d j ) r p-KV*Q ij 

where a = (p — 1) and fl = ^(p + 1). D r is the ro¬ 
tational diffusion constant and D p is a kinetic term 
also seen in prior works [31, 32, 36]. There are two 
elastic terms; De is a mean elasticity, and D$ is 
a differential elasticity, measuring the difference be¬ 
tween bend and splay energies. Finally, a fourth- 
order gradient term (with coefficient K ) is included 
in order to ensure smoothness and numerical stabil¬ 
ity. The relevant parameters for the phenomenology 
discussed are: the active force and torque (A c,r) and 
the differential elastic constant D$. In the follow¬ 
ing, we non-dimensionalize our equations by setting 
our time scale to be the rotational diffusion time, 
and our length-scale to be the diffusion length, 
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£d = y/D/D r - In all of the subsequent sections we 
will work in these dimensionless variables. 

Instabilities of the Nematic State: In the 
homogeneous limit, Eqs. (1) and (2) admit a uni¬ 
axial nematic solution with average density po > 1, 
and the order parameter Q = poSo{xx — ^1) with 

So = \J 2 ^pq+\ ~- L e t us consider spatial fluctuations 
about this state. 

Fluctuations perpendicular to the director cause 
an instability when Dq > y/2{pi - 1) ( ^^ti )- 
This instability causes phase separation into bands 
of dense ordered regions coexisting with dilute disor¬ 
dered regions when the material is near the critical 
density (p c = 1). This has been discussed in previ¬ 
ous work by us [40] and others [34-36, 38, 39]. The 
nonlinear active terms (Ac, A r and \$) do not al¬ 
ter the instability, nor do they significantly alter the 
phenomenology discussed in previous work. 

Of primary interest here is a bend instability in 
which the direction of ordering is destabilized by 
fluctuations parallel to the director. This occurs 
when a ‘bend instability parameter’ !p > 1, where 

^ = crC) ' and = 4 \/5Sy- This insta - 

bility parameter reflects a competition between the 
active torque A r and the differential elastic constant 
which is positive if the energetic cost to bend dis¬ 
tortions of the director is greater than that of splay. 
The bend instability is the primary mechanism for 
defect generation and formation of inhomogeneous 
steady states in this active system. 

In order to elucidate the consequence of this bend 
instability, we numerically explored the dynamics 
using a semi-implicit finite difference method, with 
periodic boundary conditions. Integrating from ne¬ 
matic initial conditions with small amplitude Gaus¬ 
sian noise, two states were found above the bend 
instability > 1). These were (I) an undulating ne¬ 
matic state (see Fig. 2) where the system is strongly 
ordered but the director undulates along the broken 
symmetry direction (see Fig. 1), and (II) a turbulent 
state (see Fig. 2,1) in which charge =b| disclinations 
continually form and annihilate and the +| defects 
are self-propelled as seen in [4-15]. The undulating 
nematic state was present when active convection 
Ac > 1.0 and when D$ ^ 0. In other regions of pa¬ 
rameter space the system transitioned directly into 
the defective nematic state. [41]. 

The above analysis focused on the instabilities of 
the homogeneous nematic state. Next we consider 
isotropic initial conditions and present the results for 
different values of the parameters that control the 
bend instability (A#, and po) while keeping the 
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FIG. 2. (color online) Plots showing the end state 
which forms from nematic initial conditions. Left: At 
D§ = 1.0, the homogeneous nematic state transitions 
into a defective nematic. Right: At D§ — —0.50 we find 
an undulating state at intermediate activities. 
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FIG. 3. (color online) Plots of the defect density n 
(defects per (100tA) 2 ) as a function of fj. Left: Curves 
with fixed Ac = 1.0, a range of D§ and po (scaled as 
n/(po — 1)). Defect density is greater for larger D$ when 
< 1, increases sharply when crosses 1, then satu¬ 
rates. The shaded region indicates where there was sta¬ 
tistically significant polar ordering of +| defects. Right: 
Fixed D§ = 1.0, for a range of Ac/A#, shows that when 
the active convection is comparable to active torque, de¬ 
fects vanish near the bend instability. 



other parameters fixed (Ac = D p = De = D = 1.0, 
K = 0.50, A# = 0, and Dq = 0.10) unless otherwise 
specified. 

Defect-ordered state: This nonequilibrium 
steady state occurs when defects form below the 
bend instability (vp < 1). This occurs when the 
strength of active torque is dominant over the 
strength of active convection (see Fig. 3). The prop¬ 
erties of this state are as follows: (i) Defects are 
point-like and the background is a well-ordered ne¬ 
matic (see Fig. 1). (ii) Defining the orientation of 
+ \ defects to be opposite the “comet tail’ (along the 
direction of propulsion), we find that these defects 
exhibit significant polar ordering, (iii) The degree 
of polar ordering decreases as ip increases, (iv) The 
high degree of polar ordering corresponds with long 
splay distortions which are left by defects as they 
travel (see Fig. 1). Other +| defects tend to reori¬ 
ent rather than cross this distortion-trail, leading to 
long parallel structures which are visible in the states 
with a large degree of polar ordering. 
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FIG. 4. (color online) Log-linear plot of vorticity corre¬ 
lation functions, C U (R) for fixed parameters (A r, D§ and 
po) above the bend instability (ip > 1.10). The length is 
scaled by ip, which gives a data collapse for a large range 
of parameters. 


Turbulent nematic state: The turbulent, de¬ 
fective nematic state occurs just above the bend in¬ 
stability and has the following properties: (i) Defect 
density and vorticity increase sharply for ip ^ 1. 
(ii) There is a saturation, or even a decrease, in the 
defect density (see Fig. 3) as the defects, which 
were point-like and circular near ip = 1, become 
spatially extended and the average degree of order¬ 
ing decreases, (iii) The vorticity correlation function 
C UJ (R) = (uj(0)lu(R)) / (cj 2 ) , scales with the bend in¬ 
stability parameter ip (see Fig 4), which is linear 
in the strength of the active torque A r. This dif¬ 
fers from what was found in a fluid-mediated ac¬ 
tive nematic theory [12] where vorticity scaled to 
the l/4th power of the strength of the activity, (iv) 
Assuming that the length scale for defect separation 
id scales with the vorticity (id oc A) would lead to 
a prediction that defect density scales as ip 2 . This 
seems compatible with the trend seen near the crit¬ 
ical value of the bend instability parameter, but the 
range is not large enough for a conclusive compari¬ 
son. 

Discussion: We have introduced a universal 
model of an overdamped active nematic in which 
activity enters through self-induced flows. This the¬ 
ory encompasses existing work on active nematics as 
special cases. The active torque from self-generated 
flows gives rise to a bend instability which, in turn, 
leads to defect formation. We have identified three 
nonequilibrium steady states admitted by this the¬ 
ory. The first is a defect-ordered nematic state where 
polar ordering of +1/2 disclinations emerges from 
the underlying apolar theory. The second is an un¬ 
dulating nematic state which is reminiscent of the 
“walls” of distortion in the order parameter seen 
before the onset of defective states [10, 12], or the 
distortion of the director which happens during re¬ 
laxation oscillations [42, 43] in active nematic sus¬ 


pensions. Finally we find a turbulent nematic state 
similar to that which occurs in theories of active ne¬ 
matic suspension [7-13]. 

The appearance of a defect-ordered state is of par¬ 
ticular interest, as defect-ordering was recently dis¬ 
covered in layers of active cytoskeletal filaments and 
in simulations of that system [6]. Our theory pro¬ 
vides robust predictions about when defect ordering 
will be found. The ordering occurs below the bend 
instability, i.e., ip < 1, and only for Dg > 0, which in¬ 
dicates that the energetic penalty to bend distortions 
must be greater than that of splay, and the strength 
of the active torque must be greater than that of 
the active convection. This result implies that an 
orientationally ordered fluid phase of defects may 
not occur in theories which have Galilean invariance 
(Ac = A r ). 

In [6], defect ordering in the experiment had ne¬ 
matic symmetry, while simulations displayed the 
same polar symmetry that we observe. This may 
be due to the use of periodic boundary conditions 
in both the simulations in [6] and our work, or the 
nematic ordering could be a subtle effect of fluid- 
mediated interactions in the experiment. However, 
the strong similarities between this theory and the 
experimental system, such as the sharp splay distor¬ 
tions left in the wake of traveling +^ defects leads 
us to believe that the parameter space in which this 
system operates is ip < 1, Ac < A# and Ds > 0. 
Further work is needed to substantiate this in de¬ 
tail. 

The theory that we have described contains a 
rich and diverse phenomenology which overlaps with 
other theories, simulations, and experiments of ac¬ 
tive nematics. It is generic and universal in that it 
is independent of any particular microscopic model 
and only relies on the nematic having local interac¬ 
tions. 

We thank Mike Hagan, Zvonimir Dogic, and 
Stephen DeCamp for sharing information and ideas. 
EFP and AB acknowledge support from NSF-DMR- 
1149266, the Brandeis-MRSEC through NSF DMR- 
0820492 and NSF MRSEC-1206146, NSF PHY11- 
25915 through KITP, and the HPC cluster at Bran- 
deis for computing time. EFP also acknowledges 
support through NIH-5T32EB009419 and IGERT 
DGE-1068620. 


* efputzig@brandeis.edu 
t gredner@brandeis.edu 

* baskaran@math.uci.edu 
§ aparna@brandeis.edu 










5 


[1] P.-G. De Gennes and J. Prost, The physics of liquid 
crystals , Vol. 23 (Clarendon press Oxford, 1993). 

[2] P. M. Chaikin and T. C. Lubensky, Principles of 
condensed matter physics , Vol. 1 (Cambridge Univ 
Press, 2000). 

[3] O. D. Laverntovich, Soft matter physics: an intro¬ 
duction (Springer Science V Business Media, 2003). 

[4] T. Sanchez, D. T. N. Chen, S. J. DeCamp, M. Hey- 
mann, and Z. Dogic, Nature , 24 (2013), 1301.1122. 

[5] F. C. Keber, E. Loiseau, T. Sanchez, S. J. DeCamp, 
L. Giomi, M. J. Bowick, M. C. Marchetti, Z. Dogic, 
and A. R. Bausch, Science 345, 1 (2014). 

[6] S. J. Decamp, G. S. Redner, A. Baskaran, M. F. 
Hagan, and Z. Dogic, arXiv.org , 1 (2015), 
arXiv:arXiv:1501.06228v2. 

[7] L. Giomi, M. J. Bowick, X. Ma, and M. C. 
Marchetti, Phys. Rev. Lett. 110, 228101 (2013). 

[8] S. P. Thampi, R. Golestanian, and J. M. Yeomans, 
Phys. Rev. Lett 111, 118101 (2013). 

[9] M. L. Blow, S. P. Thampi, and J. M. Yeomans, 
Phys. Rev. Lett. 113, 248303 (2014). 

[10] L. Giomi, M. J. Bowick, P. Mishra, R. Sknepnek, 
and M. C. Marchetti, Philosophical Transactions of 
the Royal Society A: Mathematical, Physical and 
Engineering Sciences 372, 20130365 (2014). 

[11] S. P. Thampi, R. Golestanian, and J. M. Yeomans, 
Physical Review E 90, 062307 (2014). 

[12] S. P. Thampi, R. Golestanian, and J. M. Yeomans, 
Philosophical Transactions of the Royal Society A: 
Mathematical, Physical and Engineering Sciences 
372, 20130366 (2014). 

[13] S. P. Thampi, R. Golestanian, and J. M. Yeomans, 
EPL (Europhysics Letters) 105, 18001 (2014). 

[14] T. Gao, R. Blackwell, M. A. Glaser, M. Better- 
ton, and M. J. Shelley, Physical review letters 114, 
048101 (2015). 

[15] X.-q. Shi and Y.-q. Ma, Nat. Commun. 4, 3013 
(2013). 

[16] D. Saintillan and M. J. Shelley, Physical Review Let¬ 
ters 99, 1 (2007). 

[17] D. Saintillan and M. J. Shelley, Physical Review Let¬ 
ters 100 , 1 (2008). 

[18] A. Doostmohammadi, M. Adamer, S. P. Thampi, 

and J. M. Yeomans, arXiv.org , 1 (2015), 

arXiv:arXiv:1505.04199vl. 

[191 V. Narayan, S. Ramaswamy, and N. Menon, Science 
317, 105 (2007). 

[20] R. Kemkemer, D. Kling, D. Kaufmann, and 
H. Gruler, Eur. Phys. J. E 1, 215 (2000). 

[21] G. Duclos, S. Garcia, H. G. Yevick, and P. Sil- 


berzan, Soft matter 10 , 2346 (2014). 

[22] R. Simha and S. Ramaswamy, Physica A: Statistical 
Mechanics and its Applications 306, 262 (2002). 

[23] S. Ramaswamy, R. A. Simha, and J. Toner, Euro- 
phys. Lett. 62, 196 (2003). 

[24] J. Toner and Y. Tu, Phys. Rev. Lett. 75, 4326 
(1995). 

[25] J. Toner and Y. Tu, Phys. Rev. E 58, 4828 (1998). 

[26] L. M. Pismen, Phys. Rev. E 88, 050502 (2013). 

[27] P. Srivastava and M. C. Marchetti, Private Comm- 
nication (2015). 

[28] P. Srivastava and M. C. Marchetti, “Instabilities 
and patterns in an active nematic film,” APS March 
Meeting (2015). 

[29] S. Mishra and S. Ramaswamy, Phys. Rev. Lett 97, 
090602 (2006). 

[30] H. Chate, F. Ginelli, and R. Montagne, Phys. Rev. 
Lett 96, 180602 (2006). 

[31] A. Baskaran and M. C. Marchetti, Physical Review 
E 77, 011920 (2008). 

[32] A. Baskaran and M. Marchetti, Phys. Rev. Lett 
101 , 268101 (2008). 

[33] S. Mishra, R. Aditi Simha, and S. Ramaswamy, J. 
Stat. Mech. Theory Exp. 2010, P02003 (2010). 

[34] F. Peruani, F. Ginelli, M. Bar, and H. Chate, J. 
Phys. Conf. Ser. 297, 012014 (2011). 

[35] A. Peshkov, S. Ngo, E. Bertin, H. Chate, and 
F. Ginelli, Phys. Rev. Lett. 109, 098101 (2012). 

[36] A. Peshkov, I. S. Aranson, E. Bertin, H. Chate, and 
F. Ginelli, Phys. Rev. Lett. 109, 268701 (2012). 

[37] A. Baskaran and M. C. Marchetti, Eur. Phys. J. E. 
35, 95 (2012). 

[38] E. Bertin, H. Chate, F. Ginelli, S. Mishra, 
A. Peshkov, and S. Ramaswamy, New J. Phys. 15, 
085032 (2013). 

[39] S. Ngo, A. Peshkov, I. S. Aranson, E. Bertin, 
F. Ginelli, and H. Chate, Physical Review Letters 
, 038302 (2013), 1312.1076. 

[40] E. Putzig and A. Baskaran, Physical Review E 90, 
042304 (2014), 1057984. 

[41] Details of the this transition, as well as the addi¬ 
tional discussion concerning the equations and lin¬ 
ear stability analysis, can be found in the supple¬ 
ment. 

[42] L. Giomi, L. Mahadevan, B. Chakraborty, and 
M. F. Hagan, Phys. Rev. Lett. 106, 218101 (2011). 

[43] L. Giomi, L. Mahadevan, B. Chakraborty, and 
M. F. Hagan, Nonlinearity , 2245 (2012). 



arXiv: 1506.03501vl [cond-mat.soft] 10Jun2015 


Supplementary Information 

(Instabilities, defects, and defect ordering in an overdamped active nematic) 


Elias Putzig, 1 ’* * * § Gabriel S. Redner, 1 ’^ Arvind Baskaran, 2, * and Aparna Baskaran 1 ^ 

1 Martin Fisher School of Physics, Brandeis University, Waltham, MA 02453, USA 
2 Department of Chemistry and Biochemistry, 

University of Maryland, College Park, MD 20742, USA 


* efputzig@brandeis.edu 

1 gredner@brandeis.edu 

* baskaran@math.uci.edu 

§ aparna@brandeis.edu 



2 


I. FREE ENERGY 


The dynamics of an equilibrium liquid crystal, which were discussed in the body of the text, can be derived 
from gradient descent on the well-known Landau-de Gennes free energy [ 1 ]: 


F = 


d d r 


- |Tr(Q 2 ) + |(Tr(Q 2 )) 2 + xQMp + 

+ 2^1 (diQki ) 2 + ~L 2 {diQik)djQjk + —LAQij(piQki)djQki 




(1) 


where the homogeneous terms give a second order phase transition and the terms with coefficients Li are 
elastic terms, associated with the energy cost of director distortions. The elastic energy cost is widely 
described in terms of splay (i^n(V • h) 2 ) and bend (A33 [n x (V x n)) ) distortion of the director. The 
parameters in ( 1 ) are related to the splay and bend coefficients as: Kn = S 2 ( 2 Li + L 2 — SX 4 ) and K 33 = 
S 2 ( 2 Li + I / 2 + 2 SL 4 ). 

The gradient-descent dynamics from the body of the paper came from taking a derivative of Eq. 1: 


= D r [a - PTrQ^Qn + 2D E V J Q, i + D,(Mi ~ - KV‘Qi 


+ * 
P 


(zQkidkdiQij + 2[dkQki\diQij — (\diQki\djQki — -Sij[dkQim ] 2 )) 


where D r = 2 / 7 , D p = — 2 7 / 7 , K = 2^/7 and the two mean elastic terms have been consolidated by using 
the tracelessness and symmetry of Q to group them in the term with coefficient De = 2 L 1/7 = 2 L 2 / 7 . The 
differential elastic coefficient, D$ = 2 L 4 / 7 , has a range — ^ and is set by the requirement that 

the equilibrium nematic should be stable to director fluctuations. Further a fourth order gradient term (with 
coefficient K ) is included in order to ensure smoothness and numerical stability. 


II. ACTIVE STRESSES AND SELF-INDUCED FLOWS 


In the body of the text we considered the dynamics of the director when the system was in an imposed 
flow,in which case the dynamical equation becomes 

(&t T ^ ’ ^)Qij — T L7~) (.Qik^kj ^ikQkj') 

L\Ejj + X(QikEkj + E ik Qkj) r 


where fand Eij are the vorticity and strain-rate tensor associated with the flow, and A is the flow-alignment 
parameter. We then postulated a form for the self-generated flow pu = — A^VQ, and the vorticity po3 = 
A*(V x VQ). This is consistent with a stress proportional to the nematic order, a = cr a Q where coefficients 
(a a , A c and A#) are positive for systems of extensile, ‘pusher’ particles, and negative for contractile, ‘puller’ 
particles. 

We construct the tensors Llij and E^ in a straightforward manner from these flows. The z- component 
of the vorticity vector gives the amount that a fluid element is rotated by the self-generated flow in two 
dimensions, and the vorticity tensor, which comes into the dynamics, is pLlij = \^ijk^k where e is the Levi- 
Civita symbol. The strain-rate tensor, which is the symmetrized gradient of the flow E^ = \{diUj — djUi) 
will play a role if the particles tend to align or tumble in shear. 

We postulated that an active liquid crystal undergoes gradient descent dynamics so as to minimize the 
energy cost of director distortions and inhomogeneities, in the local rest frame of this self-generated flow 
arising from the activity. Using the form of the flow and vorticity arising from the activity identified above, 



our dynamical equations for an active nematic take the form 
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dtQij = — 7 —I- [9kQki\diQij ——QijdkdiQki ( 2 ) 

-^(QikdkdiQji ~ QikdjdiQki - QkjdidiQki + QkjdkdiQu) 

zp 

where A s = A A c is the strength of flow-alignment, and we have simplified the flow-alignment term in the 
expression above using the symmetry and tracelessness of the nematic order tensor. 

The form of the flow-alignment term in this overdamped case is different from the form that arises in the 
overdamped limit of the the theories in which flow arises through a Stokes equation (such as [2-8]). This 
is because the phenomenological, overdamped theory which we consider here does not have incompressible 
flows. Because of this, the first order flow-alignment term in the gradient descent dynamics (A Ej-) is not 
intrinsically traceless and must be made so by hand in order to preserve the symmetric and traceless nature 
of Q, and the second order term (A [Qik^kj + EikQkj) 7 ') is nonzero, whereas it is zero in two dimensions for 
incompressible flow. The traceless version of the first order flow-alignment term is included in the functional 
derivative in Eq. 2, as it is the same as the mean elastic term: Ej- = dk(diQkj EdjQik — SijdiQki) = V 2 Qij. 


III. LINEAR STABILITY 

In the homogeneous limit, Eq. (2) and the dynamical equation for the density, d t p = DX 2 p + Dq'V'V : 
Q, admit a homogeneous uni-axial nematic state with average density po > 1, and the order parameter 

Q = poSo(xx — \ I) with So = ^ . Without loss of generality, we have picked coordinates so that the 

direction of nematic ordering is along the x-axis of our coordinate system. We parameterize the fluctuations 
about this state in the form p(r, t) = po +Sp( r, £), Q xx ( r, t) = \poSo + dQ xx ( r , t) and Q xy { r, t) = SQ xy ( r, t). 
Further, we introduce a Fourier transform X = f dre zk r A(r, t). The resulting linearized equations in Fourier 
space are 


d t Sp = —Dk 2 Sp — DqP 2 cos(2cj))5Q xx — DqP 2 sin (2 (j))5Q x 


(3a) 


d t SQx 


+ 


f Co — \D p k 2 cos (2 (f>)jSp 

^ 2 <a 0 + 2 D E k 2 + (Ds - A s)Sok 2 cos( 20 ) 

\XsSok 2 sin (2 4>)5Q xy 


Kk 4 


^5Qx 


(3b) 


d t SQ xy = - ^D p k 2 sm(2(t))5p ^ 
— ^Sok 2 sm(2(t))5Q x 


2 D E k 2 


Sok 2 cos( 2 <fi) + DsSok 2 cos( 2 <j)) + Kk 4 ^j5Q yx 


(3c) 


where 0 is the angle between the director (along the x axis in our coordinates) and the spatial gradient vector 
k. Also, we have defined parameters Co = 2 J 1 ^ ( Po ^ Q P +i 1 ) ? an d = Po — 1 fc> r notational compactness. 

The stability of these equations is best explicated by considering spatial fluctuations along (0 = 0) or 
perpendicular to (0 = |) the direction of ordering. In these two sectors, the fluctuations in the direction of 
ordering 5Q xy decouple form those of the magnitude of ordering SQ XX and the density Sp, thereby enabling 
a clear identification of mechanisms at play. 

The decoupled dynamics of SQ xy are of the following form: 


d t SQ 


xy 


= - (2D e (1 T 4>)k 2 + Kk 4 yQ xy 


(4) 


where the coefficient of k 2 is proportional to (1 — VO for fluctuations along the director (0 = §) and (1 + VO f° r 
fluctuations which are perpendicular (0 = 90°). The stability parameter, which is useful for describing the 
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dynamics when bend in unstable, is ip = , where A(po) = 4y 2 (^- 1 ) • This instability is the primary 

mechanism for defect generation and formation of inhomogeneous steady states in this active system, as was 
discussed in the body of the paper. 

Note that if the nematic was contractile, Xr is negative and hence there exists a splay instability, which 
occurs when < — 1. The analogous instability in fluid-mediated dynamics, which arises due to contractile 
stresses in the Stokes equation, has been studied in some detail [5, 7, 8]. This splay instability leads to 
defect formation, despite the fact that contractile systems enhance defect annihilation [2]. We have, however, 
focused on extensile systems for this study of the overdamped dynamics, and left the study of the contractile, 
overdamped dynamics to future work. 


IV. ADDITIONAL NUMERICAL PHENOMENOLOGY 




FIG. 1. These plots indicate whether an undulating nematic, or a turbulent, defective nematic state formed out of 
nematic initial conditions, from a point near the bend instability (if = 1.05, po = 2.0). Left Panel: A flow-tumbling 
system, with A s = —1.0. Right Panel: A flow-aligning system with A s = 1.0. The formation of defects near the 
critical density is facilitated by having a strength of active convection which is small compared to the strength of the 
active torque ( Xc/Xr < 1), and by a large energetic penalty to bend deformations, compared to the penalty to splay 
deformations (. D$ > 0). Flow-tumbling systems also facilitate the formation of defects near the critical density, but 
flow-aligning versus flow-tumbling seems to be a small effect. 
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FIG. 2. This plot shows defect density (defects per (100tA) 2 ), as a function of the bend instability parameter, 
below the bend instability (ip < 1). The shaded region shows where +| defects showed statistically significant polar 
ordering. This plot is a zoomed-in version of the plot on the left in Fig. 2 of the paper, and has parameters A c — 1-0, 
and po = 2.0. The two points indicated are systems which are shown in Fig. 1 in the body of the paper. A is the 
defect-ordered state in the upper left, and B is the defective state below the bend instability, in the lower left. 
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FIG. 3. Log-scale plots of the following approximation of the energy spectral density E(k) = ^ f d 2 Re lk ' R (v( 0) • 
v(R)). Data shows two sets of parameters. Top: (po, A r, D$)=( 2, 8,1) and Bottom: (3,10, 0.5); scaled by multiplica¬ 
tive constants for comparison. The energy density seems to scale with a power laws of k 5 ^ 3 ( k 2 ) for small wavevector, 
and k~ 5 ^ 3 (AG 2 ) and for large wavevector for the upper (lower) curve. These u-shaped energy spectra are similar 
to what is seen in mixing flow in three dimensions, and is reminiscent of that seen in bacterial turbulence and a 
reversing-rod model [9]. 




















